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Abstract 



Aims. Magnification and de-magnification due to gravitational lensing will contribute to the brightness scatter of Type la supernovae 
(SNe la). The purpose of this paper is to investigate the possibility to decrease this scatter by correcting individual SNe la using 
observations of galaxies in the foreground, without introducing any extra bias. 

Methods. We simulate a large number of SN la lines of sight populated by galaxies. For each line of sight the true magnification 
factor and an estimate thereof are calculated. The estimated magnification factor corresponds to what an observer would infer from a 
survey like SNLS. Using the simulated data we investigate the possibility to estimate the magnification of individual supernovae with 
enough precision to be able to correct their brightness for gravitational lensing with negligible bias. 

Results. Our simulations show that the bias arising from gravitational lensing corrections of individual SNe la is negligible for 
current and next generation surveys and that the scatter from lensing can be reduced by approximately a factor 2. The total scatter in 
the SN la magnitudes could be reduced by 4% for an intrinsic dispersion of 0.13 mag. For an intrinsic dispersion of 0.09 mag, which 
might not be to unrealistic for future surveys, the total scatter could be reduced by 6%. This will reduce the errors on cosmological 
parameters derived from supernova data by 4-8%. The prospect of correcting for lensing is thus very good. 

Key words, supernovae: general - gravitational lensing 



1. Introduction 



Large dedicated sup ernova surveys, such as SNLS (|Astier et al 
2006), ESSENCE (Miknai tis et all l2007t IWood- Vasey etal 



2007), and SDSS-II dFrieman et al.ll2008l) . have recently gath- 
ered data on a large number of supernovae at cosmological 
redshifts. Consequently, supernova cosmology has now reached 
a stage where systematic uncertainties, rather than lack of 
statistics, limit the ability to constrain cosmological models. 
Extinction by dust and peculiar motions of host galaxies seem 
to be the syst ematic effects wh ich currently pose the largest dif- 
ficulties (e.g. Leibundgut 2008). In this paper we investigate the 
possibility to correct for another systematic uncertainty - gravi- 
tational lensing. 

In an inhomogeneous universe like ours, light from a dis- 
tant Type la supernova (hereafter SN la) is inevitably affected by 
gravitational lensing. Weak gravitational lensing, which is what 
we will consider in this paper, refers to the phenomena which oc- 
cur when a bundle of light rays is distorted by the gravitational 
fields exerted by matter. The distortion can lead to magnification 
or de-magnification of the SN la flux. This (de)magnification 
can be described by the magnification factor, /i, which depends 
on the distribution and composition of the foreground matter. 
Both baryonic and dark matter influence the trajectories of pho- 
tons, but since the latter dominates on the scales relevant for this 
work, we only use the baryonic component as a tracer of the dark 
matter. If an observer residing in a homogeneous universe would 
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measure the flux /, an observer situated in an inhomogeneous 
universe will measure the flux / Q b s = /tt/- Due to flux conserva- 
tion, as long as we do not have multiple images, the probability 
distribution function for the magnification factor, P(n), satisfies 
the constraint (//) = J P({j,)[j,d{i = 1, i.e. the average magnifi- 
cation factor is unity. The effect of gravitational lensing conse- 
quently average out, which implies that this systematic effect can 
be controlled by observing large numbers of SNe la. In reality 
however, (//) conv erges very slowly to wards unity as the sam- 
ple size increases dHolz & Lindedl2005l) . This owes to the fact 
that the whole distribution of -P(/i), including its elusive high 
magnification tail, must be sampled, whi ch in addition to large 
numbers also requires large survey areas dCoorav et al.l l2006). 

Although the bias in principle can be beaten down by large 
numbers, gravitational lensing still contributes to the brightness 
scatter of SNe la. This contribution i ncreases with redshif t and 
is expected to be substantial at z <; 1 dHolz & Linderl2005l) . For 
precision measurements of cosmological parameters it is there- 
fore desirable to minimize this scatter by measuring and correct- 
ing for the gravitational lensing effect. 

The method we rely upon is based on modeling the dark mat- 
ter haloes, which surrounds galaxies, using the observed proper- 
ties of the galaxies. This method takes into account contributions 
to the lensing from sub-arc-minute sca les, which makes it a vi- 
able method in contrast to shear maps dDalal et alj|200l) which 
are insensitive to these scales. The possibility to correct individ- 
ual SNe la for gravitatio nal lensing using the met hod considered 
here was investigated by Gunnarsson et al. (2006), where simu- 
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lations showed that the lensing scatter could be reduced by al- 
most a factor 3 for a source at redshift z — 1.5. In this paper 
we focus on the possibility to correct for gravitational lensing 
at the lower redshifts (z < 1) accessible to the aforementioned 
large ground based surveys. The si mulations aim at m imicking 
the supernova legacy surve y (SNLS. | Astier et al.l l 2.Q06h , and fol- 
low closely the strategy by Jonss on et alj (l2008[ here after J08). 

For a correction to be useful the following two criteria must 
be met: 

1 . The correction should decrease the scatter in SN la magni- 
tudes due to gravitational lensing. 

2. Any bias introduced in the SN la magnitudes from the cor- 
rection should be negligible. 

In Sect. [2] we use simulated data to investigate the possibility 
to perform corrections satisfying the two criteria above. The 
method we use to correct for lensing is explained and applied 
to the simulated data in Sect. [3] Our results are discussed and 
summarized in Sect. [4] 

2. Simulations 

Since our simulations are very similar to the ones presented in 
more detail by J08, we only provide a brief presentation here. 

We have simulated deep galaxy catalogs with properties 
(galaxy positi ons, redshifts, lum inosities, types) based on real 
observations (Dahlen et al. 2005). For a single line of sight we 
first compute the true magnification factor, /ittrue> using the mul- 
tiple lens plane alg orithm implemented in the Q-LET package 
(Gunnarsson 2004). All galaxies are included in the calculation 
and the dark matter haloes surrounding eac h galaxy are modeled 
by truncated Navarro-Frenk-White (NFW, iNavarro et all 1 1997b 
profiles. Then we compute the magnification factor which an 
observer (with access to SNLS data) would estimate, /i ost , for 
the same line of sight. The observer would, for example, only 
see galaxies brighter than the magnitude limit (i' AB < 25.5) of 
the survey and would probably have access only to photometric 
redshifts (with a precision of oa z /n+z) — 0-03 for most galax- 
ies, but failing catastrophically for 2-6% of the galaxies). The 
content of the line of sight is therefore rather different when the 
true and estimated magnifications are calculated. Since the ob- 
server does not know what the correct model is, there is also 
a difference in the model used to compute the true and the es- 
timated magnification factor. Since the effects of gravitational 
lensing are redshift dependent, we simulate gravitational lens- 
ing of sources in the redshift range 0.2 < z < 1.1, relevant for 
e.g. the SNLS survey. 

The result of the simulations are hence pairs of true, /it ruc , 
and estimated, /j ost , magnification factors. Since most cos- 
mology fits to date have been performed using magnitudes 
rather than fluxes, we treat magnification in terms of logarith- 
mic units. Throughout the paper we denote true and estimated 
magnifications by Am truc = — 2.5 log 10 /x trU c and Am cst = 
-2.5 lo 

gio Most, respectively. 
Figure [TJ shows simulated pairs of magnification factors in 
logarithmic units. This plot shows most of the simulated pairs, 
but not the small fraction (< 0.3%) belonging to the high mag- 
nification tail with Amtruc S —0.3. The scatter in this plot 
comes mostly from the uncertainty in the model used to convert 
galaxy luminosity to mass (in our case Tully-Fisher and Faber- 
Jackson relations). Given the importance of the luminosity-to- 
mass model, we continue this investigation by considering three 
different scenarios, all shown in Fig.Q] 
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Figure 1. True magnification, Am trU c, versus estimated magni- 
fication, Awiest, for three different scenarios: no shift in halo 
masses (circles), underestimated halo masses (squares), and 
overestimated halo masses (triangles). The straight lines show 
the best fits to the no shift (solid line), underestimation (dashed 
line), and overestimation (dotted line) scenario. 



1. No shift in halo masses (circles). 

2. Underestimated halo masses (squares). 

3. Overestimated halo masses (triangles). 

If the masses of the dark matter haloes are substantially overesti- 
mated (by 50% in scenario 2) or underestimated (by 50% in sce- 
nario 3), the gravitational magnification estimates will be erro- 
neous. Our simulations aim to investigate to what extent a proper 
correction can be made for the lensing magnification given these 
uncertainties. Other sources of uncertainty, such as scatter in the 
luminosity-to-mass model, were also included for all three sce- 
narios. Clearly, the distributions of points in the plot differ for the 
three scenarios. True and estimated magnification factors are ap- 
parently correlated in all scenarios, albeit with different slopes. 
The no shift scenario has a slope (outlined by the solid line in 
Fig. QJi near unity, which reflects the fact that the model used to 
estimate the magnification is very close to the correct one. The 
underestimation (dashed line) and overestimation (dotted line) 
scenarios have slopes which differs from unity. 

From Fig. QJ the asymmetry of both the distribution of 
Am trU c and Am cst is evident. For slopes differing from unity 
this asymmetry implies an asymmetry in the errors in the esti- 
mated magnifications. The estimated magnification of a SN la 
in the high magnification tail (negative values in Fig. [TJ is likely 
to be more erroneous than a typical SN la which is slightly de- 
magnified (positive values in Fig . [TJ . 

3. Gravitational lensing corrections 

3. 1. The magnification-residual-diagram 

Since gravitationally magnified (de-magnified) SNe la should 
be brighter (fainter) than the average supernova, a correlation 
between estimated magnification and SN la brightness is ex- 
pected. According to Fig. [TJ which can be seen as an idealized 
magnification-residual-diagram, a correlation is expected even 
if the model used to estimate the magnification factors is biased. 
A real magnification-residual-diagram, where Hubble diagram 
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residuals are used instead of Amt ruc , will also be smeared by, 
e.g., measurement errors and intrinsic brightness scatter. We take 
the supernova Hubble diagram residuals to be the differences 
between observed magnitudes, m b s , and magnitudes predicted 
by a cosmological model of a homogeneous universe (we as- 
sume a flat universe, dominated by a cosmological constant, with 
Om = 0.3 and Hq — 70 km s _1 Mpc -1 ). In the following, we 
assume that the residuals have been computed using the correct 
cosmological model. The residual is then the sum of Amt rae 
and noise, Ar7i no i S0 , from, e.g., intrinsic SN la brightness scatter 
and measurement errors. We ignore other systematic sources of 
uncertainty and assume the noise to be Gaussian. 

3.2. Magnitude bias 

Because of gravitational lensing, an observer in an inhomoge- 
neous universe measures the magnitude 

TO obs = m + Am true , (1) 

rather than the magnitude to = — 2.51og 10 /, which would be 
measured if the observer resided in a homogeneous universe. 
The effect of gravitational lensing average out because the av- 
erage magnification factor is unity. However, an infinite num- 
ber of SNe la is in principle required, since the average mag- 
nification factor of a finite sample converges very slowly to 
unity as the sample size increases [see Fig. 2 in iHolz & Linderl 
d2005l) l. When magnitudes are used instead of fluxes, on the 
other hand, the effect no longer average out even in principle, 
because (log 10 /i) ^ even though (fi) = 1. Observed mag- 
nitudes are consequently biased with an amount depending on 
the redshift. The thick solid curve in Fig. [2^ shows this bias, 
("t bs) — m = (Amtrue) , whi ch is of the order 10~ 3 . 

Recently ISarkar et alj d2008l) investigated the expected bias 
in the dark energy equation of state parameter due to grav- 
itational lensing f or future SNe la surveys. In their study 
Sark ar et al] ([2008) considered both averaging ove r magnitudes 
(as we do here) and a fl ux-averaging technique dWangl l2000; 
Wang & Mukheriee 2004). In both cases they found the resulting 
bias in the equation of state parameter to be negligible, compared 
to the expected statistical uncertainty. 

3.3. Lensing corrections 

If the estimated magnification factors are used to correct for 
gravitational lensing, the corrected magnitudes would be 

m corr = TO bs - Am cst = to + e, (2) 

where the correction is characterized by the difference 

£ = ATOtrue - ATO est . (3) 

For the correction to be useful, according to the first criterion 
stated above, the dispersion in the SN la magnitudes must be 
decreased (<7 mcorr < <r mobB ), which consequently means that the 
dispersion in e must be smaller than the dispersion in Amt m e> 

i.e. Cr £ /(7Am truc < 1. 

The second criterion requires that the correction does not in- 
troduce any extra bias in the corrected magnitudes. According 
to Eq. (O this bias is given by (e) . Whether the bias is negligi- 
ble or not depends on the size of the statistical error. We require 
(e) / <TAm tIUC < 1 for a correction to be considered useful. 

Figures [2^ and [2J3 show the bias and dispersion, respec- 
tively, for e computed for the three different scenarios discussed 
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Figure 2. Bias (panel a) and dispersion (panel b) in e as a func- 
tion of redshift. The no shift scenario is represented by the thin 
solid curves. Dashed and dotted curves correspond to the under- 
and overestimation scenarios, respectively. Bias and dispersion 
of Amtrue, outlined by the thick solid curves, are also shown for 
comparison. The estimated magnifications, Am cst , were used to 
correct for gravitational lensing. 

in Sect. [2] over a range of redshifts. In this case, no noise has 
been taken into account. The scenario corresponding to no shift 
in halo masses is shown by the thin solid curves. Thin dashed 
and dotted curves outline the results for the scenarios were halo 
masses have been under- and overestimated, respectively. We 
also display with the thick solid curves the bias and dispersion 
in Am t ruc- For all three scenarios, the first criterion above is ful- 
filled, i.e. <7 e < cr Amtrue at a ^ ^shifts. At z ~ 0.8, which is the 
redshift where the SNLS distribution of SNe la peaks, the re- 
duction in the scatter is a factor 2.1, 1.9, and 1.8 for scenario 1, 
2, and 3. The absolute value of the bias, |(e)|, is < 3 x 10~ 3 
for all scenarios. If this bias is divided by erAm truc , we find 
|(e)/cAm trU cl ^ 4 x 10~ 2 , which implies successful correc- 
tions for all cases and at all redshifts according to our second 
criterion. Moreover, for all three scenarios |(e)| < |(Am t ruc)|- 
We have thereby demonstrated that such a correction for gravi- 
tational lensing would clearly decrease the scatter in the Hubble 
diagram without adding significant bias after the corrections. 

Let us now investigate the potential benefits of gravitational 
lensing corrections in the presence of realistic intrinsic bright- 
ness scatter and measurement errors. We have simulated SNLS 
like data sets, consisting of 250 SNe la each, for the three sce- 
narios. Since the results depend on the assumed errors, we have 
also varied the intrinsic brightness scatter in the simulations. To 
model the measurement errors we have use d the following fit to 
the first year SNLS data dAstier et al.ll2006l) : 

J 0.05 mag if z < 0.8 

°' crr ~ \ 0.84z 2 - 1.04z + 0.34 mag if z > 0.8. { ' 

For each simulated data set, <JAm tIUC , & e , and (e) were com- 
puted. Figure [3] shows scatter plots of (e)/<7Am t rue ( tne second 
criterion) versus <x e /<7Aro tr ue (the first criterion) for different as- 
sumed errors. As in Fig. Q] the three scenarios are indicated by 
different plotting symbols. In Fig. [3^ measurement errors, mod- 
eled by Eq. (0), and intrinsic dispersion, a- mt — 0.13 mag, corre- 
spond to the precision of first year SNLS data. A future decrease 
in the intrinsic dispersion is not to unrealistic. We therefore show 
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in Fig.[3p results obtained for the same measurement errors, but 
with (Ti n t = 0.09 mag. FigurefS}; shows results for <Tj llt = 0.09 
mag and negligible measurement errors, i.e. <7 Grr = mag. For 
comparison we also show results in Fig. [3jl for an intrinsic dis- 
persion of only 0.05 mag and er crr = 0. The distributions of 
points in the scatter plots are rather similar irrespective of the 
scenario. A value of cr e /<7Am truo less than unity, i.e. to the left 
of the vertical lines, indicates a successful correction according 
to the first criterion. The number of simulated data sets failing 
to meet the first criterion, i.e. points to the right of the vertical 
line, decreases as the errors decrease. In Fig. [3^ the correction 
fails for 4%, 1%, and 10% of the synthetic data sets for the no 
shift, underestimation, and overestimation scenario. These num- 
bers drop to less than one percent for Fig.JSJi. 

From Fig.[3]it is also evident that the bias after the correction 
has been performed is small; for all cases | (e) / <7Am tluc I < 0.25. 
According to Fig. [3] corrections for gravitational lensing are 
likely to reduce the scatter in SN la brightness without intro- 
ducing any harmful bias. 

Figure 0] shows a projection of Fig. [3] focusing on the dis- 
tribution of the ratio a e j 'o"Am trU c- Figures 0^, 0b, and [4]: show 
results obtained for the no shift, under-, and overestimation sce- 
narios. Solid and dashed curves correspond to measurement er- 
rors modeled by Eq. (0]l and intrinsic dispersion of 0. 13 and 0.09 
mag, respectively. Dotted and dash-dotted curves correspond to 
negligible measurement errors (cr err = 0) and intrinsic disper- 
sion 0.09 and 0.05 mag, respectively. In TableQ] the character- 
istics of the distributions shown in Fig.0]are summarized. From 
Fig- IH it is clear that the distributions obtained with systematic 
shifts in the halo masses (scenario 2 and 3) are rather similar 
to the distributions corresponding to the no shift scenario. From 
Table [T] we see that 0.96 < (cr e /oAm truc ) < 0.97 when the 
intrinsic scatter is <7j n t = 0.13 mag, which implies that the av- 
erage effect of correcting for gravitational lensing is to decrease 
the total scatter by ?>-\%. A realistic future intrinsic dispersion 
of (Ti n t = 0.09 mag would allow a decrease in the total scatter by 
5-6%, for measurement errors modeled by Eq. and by 7-8% 
if measurement errors could be reduced to a negligible level. 

We have also investigated the possible improvements on cos- 
mological parameter estimation which could be gained from 
gravitational lensing corrections. Using the Fisher information 
matrix, we computed confidence ellipses in the (fi^i, -plane 
assuming a flat universe. The parameters 17m and w denote the 
dimensionless matter density of the universe and a constant dark 
energy equation of state parameter, respectively. In this analy- 
sis, a low and a high redshift SN la data set were used. The 
high redshift data set consists of 250 SNe la with redshif ts and 
meas urements errors similar to first year SNLS data (Astie retail 
2006). The purpose of the low redshift data set, consisting of 44 
SNe la unaffected by lensing, is to anchor the Hubble diagram. 
Redshifts and uncertainties we re similar to the low redshift data 
set used in lAstier et alj ([2006). Errors due to gravitational lens- 
ing before and after correction were modeled using the curves in 
Fig. [2] Our treatment of the errors is hence simplified, because 
these curves correspond to the dispersion and does not take the 
asymmetry of the distributions into account. Table [2] shows the 
relative improvement in the area of the confidence level ellipses 
due to gravitational lensing corrections, ^4 C orr /^tensed- The scat- 
ter due to gravitational lensing is rather small compared to intrin- 
sic brightness dispersion and measurement errors. Nevertheless, 
corrections for gravitational lensing would result in improve- 
ments of 4-8% for realistic values of the intrinsic dispersion and 
~ 30% for a very optimistic scenario where other sources of 
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Figure 3. Scatter plot of (e) /<JAm tIUO versus cr e /oAm truc f° r 
simulated SNLS like data sets consisting of 250 SNe la each. 
Circles, squares, and triangles represent the no shift, under-, and 
overestimation scenarios, respectively. The different panels cor- 
respond to different intrinsic brightness dispersion and measure- 
ment errors: Panel a) <jj nt = 0.13 mag, tr orr given by Eq. (0); 
Panel b) <7j nt = 0.09 mag, cr crr given by Eq. (0); Panel c) 
CTint = 0.09 mag, a CII = mag; Panel d) a- m t = 0.05 mag, 

a^r = mag. 



error are negligible and the supernovae are calibrated to 5% ac- 
curacy. 
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Figure 4. Probability distribution functions of CT e /crAm tluc f° r 
different scenarios and errors. The distributions were obtained 
for simulated SNLS like data sets consisting of 250 SNe la. 
Panel a, b, and c correspond to the no shift, under-, and over- 
estimation scenario. The curves show results for different as- 
sumed errors: Solid curves correspond to er crr given by Eq. (0| 
and (Ti n t = 0.13 mag; Dashed curves correspond to <j err given 
by Eq. (0} and a mt = 0.09 mag; Dotted curves correspond to 
0"err = and (Tint = 0.09 mag; Dash-dotted curves correspond 
to (Ten = and (T; nt = 0.05 mag. 
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Table 1. Characteristics (mean and root mean square) of the dis- 
tribution of <T e /CTAm truc for me no shift, under- and overestima- 
tion scenarios and different assumed errors. The distributions 
were obtained from 10000 simulations of SNLS like data sets 
consisting of 250 SNe la each. 



(Tint cr err Scenario 1^ Scenario 2^ Scenario 3^ 

(mag) (mag) Mean RMS Mean RMS Mean RMS 

013 EqU 096 O02 097 O02 097 OOF 

0.09 Eq. 0) 094 0.03 0.95 0.02 0.95 0.04 

0.09 0.92 0.04 0.93 0.03 0.93 0.04 

0.05 0.80 0.06 0.83 0.05 0.83 0.07 



" No shift in halo masses. 

* All halo masses are underestimated by 50%. 

c All halo masses are overestimated by 50%. 

Table 2. Relative improvement in confidence level contours in 
the (17m, w)-plane from corrections for gravitational lensing. 
The results presented in this table were obtained using a Fisher 
matrix analysis for a SNLS like data set consisting of 250 high 
redshift SNe la and 44 low redshift SNe la (assumed not to be 
affected by lensing). 
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Figure 5. Correction coefficient B as a function of redshift for 
three simulated scenarios. For the first scenario (solid curve) 
there is no systematic shift in halo masses. For the second 
(dashed curve) and third (dotted curve) scenario halo masses 
were systematically under- and overestimated by 50%. 



o"int Terr Scenario 1^ Scenario 2_ Scenario 3^_ 

(mag) (mag) ^4corr/^llcnscd ^4corr/^4lcnscd ^.corr /^llcnscd 

0T3 EcpH 095 091 096 

0.09 Eq. © 0.92 0.93 0.93 

0.09 0.84 0.87 0.87 

0.05 0.69 0.74 0.74 



" No shift in halo masses. 

b All halo masses are underestimated by 50%. 

c All halo masses are overestimated by 50%. 

3.4. Correction coefficient 

In the next step, we could try to improve the correction by con- 
sidering a model which incorporates the different slopes exhib- 
ited in Fig.[T] Since the points in Fig.[T]appears to cluster around 
straight lines with different slopes, a linear correction in ATO es t, 

Amf st = BAm cst , (5) 

is the simplest approach. The £?-coefficient is related to how bad 
we are at estimating the luminosity-to-mass relation, and cor- 
recting with B is the most straightforward way to remedy this. 
In J08 it was shown that since B is sensitive to the normalization 
of halo masses, we can constrain galaxy masses using the fitted 
value of B. 

Since the effects of gravitational lensing are redshift de- 
pendent, the correction coefficient could change with redshift. 
Figure[5]shows the best fit values of B to a large number of sim- 
ulated pairs of ATO cst and Ara trac for our three scenarios as a 
function of redshift. An optimal correction would, of course, use 
values of B obtained at different redshifts, but would be difficult 
to derive for a limited data set. Fortunately, the correction coef- 
ficient B stays rather constant with redshift, which implies that 
a slope fitted to a data set spanning a range of redshifts should 
be useful. 

Figure |6l which is is analogous to Fig.|2j shows dispersion 
and bias as a function of redshift after corrections have been per- 
formed using the best fitting correction coefficient, B, to the en- 
tire data set, which spans a range of redshifts. Equation (0 can 
not remove the bias at all redshifts, but can bring the bias for the 



under- and overestimation scenario down by ~ 25% and ~ 70%, 
respectively. 

Including the correction coefficient in the corrections have 
a negligible impact on (e) /erAm truc f° r our simulated data sets 
(see Sect. 13.3b . but make a difference to the distributions of 
Oe/fAmcsf Figure [7] which is similar to Fig. |4] shows the 
distributions for the no shift scenario when corrections have 
been performed using Eq. (0. Only the distributions corre- 
sponding to scenario 1 is shown, because the distributions be- 
longing to the three different scenarios are indistinguishable. 
Equation (0 hence eliminates the effect of the biased halo 
masses. Furthermore, all distributions are pushed to the left of 
the vertical line (er e /<7Am cst = 1)- The inclusion of the correc- 
tion coefficient B, thus reduces the corrected gravitational lens- 
ing dispersion down to the level of scenario 1, even if the halo 
masses are originally over- or underestimated. 

It also improves cosmological constraints in the (Qm,w)- 
plane to the same level as for the no shift scenario without cor- 
rection coefficient (see Table 0. No further improvement in the 
cosmological constraints can be achieved for the no shift sce- 
nario by Eq. (0. 

3.5. Judging the correction 

Correcting for gravitational lensing can of course only be jus- 
tified if there is a correlation between residuals and estimated 
magnification. The presence or absence of a correlation after the 
correction has been performed can be used to judge the success 
of the correction. If the correction was successful, there should 
no longer be any correlation. On the other hand, if the correction 
was not justified we might have introduced a spurious correla- 
tion. 

Although the slope differ for the different scenarios plotted 
in Fig.Q] the linear correlation coefficient does not. The correla- 
tion coefficient, r, depends on the scatter, which is the same for 
all three scenarios. We have simulated a large number of data 
sets and computed the correlation coefficient before and after 
the correction. Before the correction we consider the correlation 
between Am cst and Hubble diagram residuals. After the correc- 
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Figure 6. Bias (panel a) and dispersion (panel b) in e as a func- 
tion of redshift. The no shift scenario is represented by the thin 
solid curves. Dashed and dotted curves correspond to the under- 
and overestimation scenarios, respectively. Bias and dispersion 
of Amtruei outlined by the thick solid curves, are also shown 
for comparison. To correct for gravitational lensing the formula 
Amf st = BAm cst with B corresponding to the value outlined 
by the straight lines in Fig. Q] was used. 
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Figure 7. Probability distribution functions of a £ / 'a A mtlue for 
the no shift scenario and different errors. The distributions were 
obtained for simulated SNLS like data sets consisting of 250 
SNe la corrected for gravitational lensing using the formula 
Amf st = BAm cst . The curves show results for different as- 
sumed errors: Solid curves correspond to a eir given by Eq. 
and (Ti n t = 0.13 mag; Dashed curves correspond to a crr given 
by Eq. (01 and er int = 0.09 mag; Dotted curves correspond to 
Corr = and <7i nt = 0.09 mag; Dash-dotted curves correspond 

tO (Terr = an d °mt = 0.05 mag. 

tion the correlation between Am cst and e = Am tmc — Am es t 
is considered. For all three scenarios, the initial correlation co- 
efficient is (r) ~ 0.27. The uncertainty in the correlation co- 
efficient decreases as <r r ~ 1.3/ yN as the number of SNe la 
increases. After applying the correction, we find (r corr ) to be ap- 
proximately —0.02, 0.09, and —0.10 for the no shift, under-, and 
overestimation scenario. For the under- and overestimation sce- 



narios, where the magnifications are not correctly estimated, the 
correction hence results in weak spurious correlations. Whether 
r and r corr are significantly different or not depends on the sam- 
ple size. For the no shift scenario, the difference would be signif- 
icant (above 3er) for N <; 200. If the correction coefficient was 
taken into account, Am cst and e = ATO trU c — Am es t would be 
completely uncorrected since B was fitted to the data, making 
this test inappropriate. 

4. Discussion and summary 

Future SN la data sets are anticipated to be large, which could 
justify flux averaging as a method to overcome gravitational 
lensing bias. However, such large homogeneous data sets will 
probably allow for a reduction of systematic uncertainties and 
the use of subsets of SNe la - or the discovery of new corre- 
lations with peak brightness - may lead to large reductions of 
the intrinsic brightness dispersion. Under such circumstances, 
gravitational lensing corrections may play an important role for 
precision cosmology. 

In this paper we have investigated the possibility to perform 
corrections for gravitational lensing of individual SNe la that re- 
duces the scatter with negligible bias. Since one of the largest 
uncertainties in the estimation of magnification factors is the 
relation between luminosity and mass of galaxy haloes in the 
foreground, we have studied three different scenarios: no shift 
in halo masses, underestimation by 50% for all halo masses, 
and overestimation by 50% for all halo masses. Our simulations 
show that for all three scenarios, the scatter due to gravitational 
lensing can be reduced by roughly a factor 2 for a SNLS like data 
set. Also, the correction will be useful even if we do not know 
the luminosity-to-mass relation for the lensing galaxies to bet- 
ter than 50%. Any apprehension that gravitational lensing cor- 
rections would lead to bias, thus appears to be unfounded. For 
simulated SNLS like data sets we found |(e)/erAm trU el ^ 0.25. 
The bias (e) is also expected to vary with redshift, but at a much 
smaller level ( | (e)/cr A m truc | < 4 x 10~ 2 ). 

As expected, the no shift scenario decreases the scatter more 
than the under- and overestimation scenarios. However, a sim- 
ple correction coefficient, B, parameterizing the slope in the 
magnification-residual-diagram, can increase the performance 
of the corrections of the under- an overestimation scenarios 
to the same level as the no shift scenario. Including the B- 
coefficient in the correction can hence eliminate the effect of 
bias in halo masses. 

Correcting for gravitational lensing could reduce the total 
scatter in the SN la magnitudes with 4% for an intrinsic disper- 
sion of (Ti nt =0-13 mag and errors similar to the ones obtained 
for SNLS dAstier et al.l2006l) . For cr int = 0.09 mag, which might 
not be to unrealistic for future surveys, the total scatter could be 
reduced by 6%. Using Fisher matrix analysis we find that the 
size of the error ellipses in the (J7m, w)-plane can be reduced by 
4-8% for realistic measurement errors and realistic values of the 
intrinsic brightness scatter. 

Gravitational magnification of SNe la is rather independent 
of cosmological parameters. The cosmology dependence on the 
corrections should thus be rather small. However, when a more 
elaborate correction, such as Eq. (0, is employed which requires 
a fit to data we have to be more careful. The residuals used 
depends on the cosmology and a proper correction must take 
this into account. One solution would be to iteratively compute 
the slope in the magnification-residual-diagram and estimate the 
cosmological parameters. In passing, we note that the gravita- 
tional lensing magnification we consider here affects all points 
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of the SN la light-curve by the same amount, making it possible 
to apply the corrections after the light-curve fit. 

Corrections for gravitational lensing could also be impor- 
tant for other distance indicators. Gravitational waves emitted by 
chirping binary systems could - if detected - be used to obtain 
very accurate luminosity distances dSchutz| [l986). Furthermore, 
if the redshift of the binary could be measured from an optical 
counterpart, these binary system could be standard sirens, the 
gravitational wave analogs of standard candles. These standard 
sirens are unaffected by most systematic uncertainties which 
plague standard candles and they might provide distances with 
relative accuracy < 1%. Standard sirens and standard can- 
dles both suffer extra dispersion due to gravitational lensing. 
For standard sirens this effect is much more important than 
for standard candles, since the other uncertainties are so small. 
Gravitational lensing will thus degr ade the power of chirpi ng bi- 
nary systems as distance indicators jHolz & Hug hes 2005). The 
potential improvem ent of correcting sta ndard sirens for lensing 
was investigated in Jon sson et al.l (120071) . Corrections could re- 
store some of their power. The results found here that gravita- 
tional lensing corrections can be unbiased could thus be of im- 
portance also for future precision gravitational wave cosmology. 
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